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Abstract 

We consider the process of scattering of Galactic cosmic-ray electrons and protons off of dark 
matter with the radiation of a final-state photon. This process provides a novel way to search for 
Galactic dark matter with gamma rays. We argue that for a generic weakly interacting massive 
particle, barring effects such as co-annihilation or a velocity-dependent cross section, the gamma- 
ray emission from cosmic-ray scattering off of dark matter is typically smaller than that from dark 
matter pair-annihilation. However, if dark matter particles cannot pair-annihilate, as is the case 
for example in asymmetric dark matter scenarios, cosmic-ray scattering with final state photon 
emission provides a unique window to detect a signal from dark matter with gamma rays. We 
estimate the expected flux level and its spectral features for a generic supersymmetric setup, and 
we also discuss dipolar and luminous dark matter. We show that in some cases the gamma-ray 
emission might be large enough to be detectable with the Fermi Large Area Telescope. 
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I. INTRODUCTION 



Indirect dark matter detection is the search for the annihilation or decay products of 
particle dark matter — instead of directly detecting the dark matter, indirect detection 
aims at pinpointing the non-standard origin of a cosmic-ray or a gamma-ray signal based 
upon the messenger's spectral or angular features (see the relevant chapters of Ref. p] for a 
collection of recent reviews). Examples of indirect detection with gamma rays are searches 
for a peculiar gamma-ray signal from nearby dwarf spheroidal galaxies (see e.g. [2]), from 
nearby clusters (e.g. [3H5]), from the center of the Galaxy (e.g. [6]), in the extragalactic 
gamma-ray background (e.g. [HE]) or in the diffuse background angular distribution (e.g. 
[9]); high-energy neutrinos from the center or of the Sun or of the Earth (see e.g. [TQl [11]) , the 
search for features in the spectrum of positrons or antiprotons (e.g. |12j). and the search for 
antidetuerons (e.g. [321 [33]) are further examples of indirect dark matter detection channels. 

This is an especially exciting point in time for the field of indirect dark matter detection. 
The Fermi Large Area Telescope (LAT) is well into its third, successful year of science 
operations. Crucial aspects of the gamma-ray sky for indirect dark matter searches, such as 
the issue of the galactic and extra-galactic diffuse background, are undergoing close scrutiny 
with significant progress well under- way [T^l HE]- Research and development is also in 
progress for the next generation gamma-ray detector, a ground-based Cherenkov Telescope 
Array (CTA), whose significantly lower low-energy threshold will enormously benefit its use 
for dark matter searches. The IceCube neutrino telescope has recently been completed, 
with a thickly-instrumented portion (DeepCore) that makes it an even better dark matter 
detection instrument than originally designed [11]. Finally, the AMS-02 payload has recently 
been successfully deployed on-board the International Space Station, and science results will 
soon be delivered. 

The weakly-interacting massive particle (WIMP) paradigm for the nature of dark matter 
is a highly compelling scenario, as it relies on generic weak-scale physics and on the well- 
known thermal freeze-out mechanism for the production of a relic with just the right density 
to match the inferred dark matter abundance in the universe. However, the so-called WIMP 
miracle does not provide an explanation to why the universal density of baryons Qb (in units 
of the critical density, fie = Pb/Pcth) is tantalizingly close to the density of dark matter DM , 
so much as to pose the veritable "coincidence problem" of explaining why Qdm/^b ~ 5. 
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A general class of models posits that the dark matter relic density is in fact set by an 
asymmetry between dark matter and anti-dark matter particles, in full similarity to the 
baryonic relic density (e.g. [TTH20] ). If the two asymmetries are related to each other (for 
example if there exists a "transfer mechanism" between the dark matter and the baryonic 
sector), then generically m DM ps c (f2 DM /f2 B ) m P ~ 5 GeV, with m p the proton mass and c 
a model-dependent number of order unity. These models have been dubbed "asymmetric 
dark matter" models in Ref. [21], and received a considerable amount of attention recently, 
in part because of direct detection results that might point to a light (with a mass of a few 
GeV) dark matter candidate. 

While asymmetric dark matter models give rise to an interesting phenomenology (see 
e.g. [2E 122]), including new contributions to proton decay, and, as mentioned, can even 
explain anomalies in direct dark matter experiments (e.g. [23]), by definition, asymmetric 
dark matter does not annihilate (or does so at a very suppressed rate, should there be a 
left-over symmetric component). One might thus naturally wonder whether indirect dark 
matter detection is or not entirely doomed in asymmetric dark matter scenarios (see however 

[2H2U). 

In this paper we consider a process that was, in an entirely different astrophysical setting, 
suggested some time ago by Bloom and Wells [28J. There, the Authors suggested that high- 
energy electrons in the jets of active galactic nuclei (AGN) could up-scatter dark matter 
particles, possibly radiating a photon in the final state. Estimating, with the techniques of 
effective field theory, the relevant cross section, and estimating the flux of electrons and the 
density of dark matter particles, the authors concluded that (1) the level of the signal was 
several orders of magnitude below the sensitivity of the LAT, and (2) no specific "smoking- 
gun" spectral features were present in the signal. 

In Ref. [29] we recently studied the same process originally proposed in [28], improving 
on the previous analysis in a number of ways. First, we carried out the complete elementary 
process calculation, including resonance and threshold effects, for two particle dark matter 
models (supersymmetric neutralinos and the lightest Kaluza- Klein of universal extra dimen- 
sions) for both cosmic-ray electrons and protons in the AGN jet; Second, we re-evaluated 
the AGN jet components and the dark matter density in view of recent observational and 
theoretical- modeling results. We concluded that the process can be resonant, and that the 
resonance entails on the one hand a larger final-state photon flux, and on the other hand 
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a distinctive spectral feature, corresponding to a certain, well-understood kinematic condi- 
tion. As a result, we predicted that for some dark matter models the process of dark matter 
scattering off of AGN jet particles can be observable by the Fermi LAT. 

Here, we observe that our own Galaxy is populated by cosmic rays (hadrons and heavy 
nuclei, electrons, positrons and neutrinos) that might interact with the Galactic dark matter 
- a process that might be "visible" with a gamma-ray telescope such as Fermi LAT if it 
radiates a photon in the final state. Evidently, the elementary processes we considered in 
detail in Ref. [29] are precisely the same for the case of scattering off of Galactic cosmic 
rays (albeit a few technical differences will be of importance, as we point out later). A 
further possibility, that we leave for future work, is that dark matter-cosmic ray scattering 
lead to inelastic processes producing e.g. neutral pions, eventually resulting in gamma rays. 
Note that inelastic interactions between cosmic-ray protons and dark matter particles have 
been also recently considered in Ref. [30J and [3T] as a possible explanation to the origin 
of the observed "knee" in the spectrum of galactic cosmic rays, for dark matter models 
in the context of extra dimensions with a low fundamental scale of gravity. Limits on 
the interaction cross section of dark matter off of hadronic cosmic rays have also been 
considered in Ref. |32j in the context of models of strongly interacting dark matter; the 
strongest constraints were found to arise from effects on the predicted abundances of light 
elements in big bang nucleosynthesis. 

In this study we first consider, in Sec. [TTJ, the process of dark matter scattering off of cosmic 
rays in the Galaxy for a generic WIMP model, and we compare the expected luminosity 
from final-state radiated photons with what expected for WIMP pair-annihilation; we also 
consider the same process for two additional dark matter models that received some attention 
recently: dipolar dark matter and "luminous" dark matter. We then present the core of our 



numerical results in Sec. Ill, and separately give details for the case of cosmic-ray electrons 



(III A) and protons (III B). We present a discussion and our conclusions in the final Sec. IV 



II. DARK MATTER PAIR ANNIHILATION VERSUS COSMIC RAY-DARK 
MATTER SCATTERING 

In this Section we estimate, for generic WIMP models as well as for dipolar and luminous 
dark matter models, the expected luminosity for the process where a photon is emitted in 
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the final state of a cosmic-ray scattering off of dark matter. 

We start by comparing, for the case of (symmetric) WIMP dark matter models where 
pair-annihilation can occur, the luminosity of the final-state photon emission from cosmic 
ray scattering off of dark matter to the well-known dark matter pair-annihilation rate. Both 
processes can be thought of as occurring on the same target particles: Galactic dark matter. 
We thus construct an order-of-magnitude estimate for the relevant rates, per unit target. 
These are defined here as follows: 

• for dark matter annihilation, 

IxX = a XX n X U rel) 

with a xx the pair annihilation cross section, n x the number density, and v Te \ the relative 
dark matter particle velocity in the Galaxy 

• for cosmic ray-dark matter scattering, 

A4>e 



dE 



1 GeV 



with the locally measured differential cosmic ray electron flux, averaged over angles, 
and cr ex -> eX7 the relevant scattering cross section. 

We estimate a weak-scale pair-annihilation cross section 

2 

a xx ~ — — ~ (few) x 1(T 36 cm 2 , 



m EW 



and we employ a local dark matter density p x ~ 0.3 GeV cm -3 and a dark matter mass 
of 100 GeV to estimate n x ~ 3 x 10~ 3 cm -3 . Finally, we take v Te \ ~ 10~ 3 c. With these 
numbers, we find 

q xx ~ (few) x lO" 31 - 

For the case of cosmic ray-dark matter scattering we employ 

1 



dE 



~ 4 x 10" 1 - 



1 GeV Cm 2 s' 



as measured and reported by the LAT collaboration [33] (we consider cosmic-ray electron 
energies around 1 GeV given that we are interested in final-state gamma rays with energies 
in the Fermi LAT range, thus above 0.1 GeV, and because at larger energies the cosmic-ray 
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spectrum declines very steeply; we thus expect the energy range that contributes the most 
to the signal we are interested in to be in the GeV range), and we estimate 

a 3 

°ex^ex7 ~ 2 ~ (few) X 10 38 CHI 2 , 

m EW 

with possible resonances enhancing these numbers by 2-4 orders of magnitude, according to 
the findings of (29] (see e.g. Fig. 3). As a result, 

q ex ^exi ~ 10~ 38 - -> (few) x 1(T 34 -. 

s s 

This calculation only employs the local dark matter density and cosmic-ray electron flux: 
both these quantities are larger the closer one gets to the center of the Galaxy. This order of 
magnitude estimate, however, implies that in the generic case where the dark matter pair- 
annihilates and there are no symmetries or kinematic effects suppressing the natural pair- 
annihilation rate for a weak-scale weakly-interacting particle, the effect under investigation 
here is suppressed by 3 to 7 orders of magnitude compared to pair-annihilation. Notice also 
that the number of photons produced in each annihilation event is typically on the order 
of 10-100, as opposed to the single photon that scattering off of cosmic ray would produce 
per event (an exception to this is given by inelastic processes leading to neutral pions in the 
final state). Since the detectable quantity here is not the event-rate, but the actual gamma 
rays produced in the interactions, the suppression in the predicted gamma-ray flux is thus 
generically larger by a further 1 to 2 orders of magnitude. 

If, however, effects such as co-annihilation or velocity dependence of the pair-annihilation 
are in place, the two processes could produce a comparable flux of gamma rays. In addition, 
if dark matter cannot pair-annihilate because of discrete symmetries or because of lack of 
annihilation partners, such as in asymmetric dark matter models [21] , then scattering off of 
cosmic ray electrons with a final-state photon emission might be one, if not the only way to 
detect a signal from dark matter with gamma-ray telescopes. 

It is interesting to investigate two more dark matter scenarios that have received some 
attention recently: dipolar dark matter (see e.g. |34j). and luminous dark matter [35]. The 
former arises as dark baryons in hidden-sector gauge mediation supersymmetry breaking 
(see e.g. [36]), or in technicolor models (see e.g. [37T440] ). While the relic density of this 
class of dark matter candidates is produced non-thermally, often a mechanism relates the 
dark baryon to the Standard-Model baryon density, thereby providing an explanation to 
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the coincidence problem S^dm ~ [IS]- Luminous dark matter, on the other hand, is an 
interesting proposal in light of recent direct detection results [35J. 

We estimate the dark matter scattering off of cosmic ray electrons for the case of dipolar 
dark matter assuming that the dark matter is a fermion ip of mass tbdm ~ 100 GeV, with 
a dipolar interaction cross section of the type 

g _ 

Aiipolc = l^^^F^v 

We estimate the cross section off of cosmic-ray electrons with a typical energy scale around 
1 GeV, and consider for simplicity the internal bremsstrahlung approximation of Eq. (17) of 
Ref. |29j for the final state photon radiation, with an outgoing photon energy on the GeV 
range. Notice that this process is never resonant. We obtain the estimate 

vy 

A 



This entails that unless the dark sector scale A is very low and the dark matter is very light 
the dipolar dark matter case is less promising than the WIMP case. 

In the context of luminous dark matter [35], the dipolar operator connects two Majorana 
fermions, Xg (t ne dark matter particle) and Xe, with m Xg ~ 1 GeV and m Xe = m Xg + S and 
5 ~ 3.3 keV chosen to provide an explanation to the observed DAMA/LIBRA spectrum. 
The relevant interaction term is 

Awn = JTXgV^XeF^ + ll.C. 

Here the picture is one where scattering of Xg °ff °f cosmic ray via the dipolar interaction 
operator produces meta-stable Xe particles, subsequently decaying into Xe — > Xg + 7 with 
the photon energy in the final state E 1 ~ 5 ~ 3.3 keV in the X-ray regime. 

We estimate the cross section for the process ex g — > eXe for E e ~ 1 GeV to be on the 
order of 



(1 rp t r\ 2 
) cm 2 . 
A / 



This cross section would produce a rate of 3.3 keV X-ray photons per unit dark matter 
target of 

-36 



9e X9 ^eX 9 7 ~ 10 I J 



1 TeV\ 2 1 
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Considering a Galaxy similar to the Milky Way, and taking a total mass of Mmw ~ 1O 12 M , 
this leads to an X-ray luminosity of about 

M MW in 33 1 m 

" l Xg b 

which appears to be too small to be detectable if compared, for example, to the typical 
unresolved X-ray emission for a Milky Way-type galaxy of 

pT'CT 1 

Lx-ray ~ (few) x 10 39 — ~ 10 47 - at 1 keV. 

s s 

Notice that the effective target number density is likely over-estimated in Eq. 0, as one 
should really only consider dark matter particles within the diffusion region where Galactic 
cosmic rays are confined. 

One can also compare our estimate in Eq. with the luminosity that would be produced 
by sterile neutrinos with a lifetime r ~ 10 31 s and a mass of a few keV, which lie at the verge 
of the sensitivity of current X-ray data (see e.g. [H]). This case is relevant here since sterile 
neutrinos at that mass would produce exactly the same X-ray feature, a monochromatic line 
at an energy of a few keV. We estimate that the corresponding Milky Way X-ray luminosity 
would be on the order of 

£lim _ 1Q 44 I 
S 

We thus infer that the luminosity we estimate for luminous dark matter excitation from 
cosmic-ray scattering, Eq. ([!]), is well below the current X-ray telescope sensitivity. 

In what follows, we shall focus on the case of WIMPs, motivated by the estimates above, 
and by the fact that if resonances are present, those can produce interesting and potentially 
revealing spectral features, such as those we found in the case of scattering of dark matter 
off of AGN jet cosmic rays 



III. COSMIC RAY - WIMP DARK MATTER SCATTERING 

In this section, we use the analytical results of Ref. on the neutralino-electron and 
neutralino-proton cross section with the radiation of a photon in the final state to calcu- 
late the diffuse gamma-ray galactic emission from dark matter scattering off of cosmic-ray 



electrons (sec. Ill A) and off of cosmic-ray protons (sec. Ill B ) 
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A. Electrons 



We want to compute the spectrum of photons coming from the scattering of cosmic rays 
off of the dark matter in the Milky Way. We shall consider the electrons first. The photon 
flux is given by 

dN - 1 f f dxj> d 2 a 

where 



dE^ ^'^ M X J "J a dE e dQ^dE 7 



J = M / Mf& * ej ( e )> ( 3 ) 



r 2r e i 

m = / ds p(r(s,9))f(r(s,e)), 

Jo r eP& 

and /(r) = e~ r//r ° /e~ r ®' r ° (see e.g. Ref. [12] and references therein) is included to take into 
account the fact that the cosmic ray flux is larger in the vicinity of the galactic center. In 
the expressions above, f2 7 is the angle between the emitted photon and the incoming cosmic 
ray. To account for cosmic rays coming from every direction we integrate over the full solid 
angle f2 7 . On the other hand, Afl is the angle of observation from earth. We choose to 
focus on the direction towards the Galactic Center (GC), where the dark matter density 
and the cosmic-ray density are highest; given the typical angular resolution, for the energies 
of interest, of the LAT [13], we choose an angular region of 9 = 1° around the GC, which 
corresponds to AQ ~ 10~ 3 . As far as the dark matter distribution p(r) is concerned, we use 
here the spherically-symmetric Einasto profile |Sj 



p(r) = p s exp 
The cosmic ray electron flux is taken to be 



r_2 




1 


[ a 







(4) 



# • (5) 



dE e \GeV 

for E e between ~ 0.1 and ~ 1000 GeV. The numerical values for the various parameters we 
employ in the numerical evaluations of the Equations above are summarized in Table [Tj 

The last ingredient is the differential cross section d ^f^ E , in order to compute which, one 
has to specify a particle dark matter model. For the sake of illustration, we consider here 
the Minimal Supersymmetric Standard Model (MSSM), with the lightest neutralino as the 
dark matter candidate. 
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FIG. 1: Schematic diagram for electron (quark) - neutralino scattering proceeding through an s- 
channel. 

1. Case study: Neutralino dark matter 

We recently carried out the complete calculation of the differential cross section in 
Ref. [2H], using the collinear approximation, i.e. taking the photon to be collinear with 
the final electron. We relax here, however, that assumption, while still only considering 
diagrams where the selectron is exchanged in the s-channel (see Fig. [TJ, since those are the 
only possibly resonant channels. Using the same notation as in [29J, it is useful to remind 
the reader about the values of the coupling between electron (quark) - selectron (squark) - 
neutralino: a L = V2qg tanOw and or = l/2a^, with q the electric charge in units of e, g 
the Standard Model SU{2) gauge coupling and 6 W the Weinberg angle. R and L refer to 
the handedness of the electron or quark. 

In Ref. [29], the directions of the initial electron and the final photon were fixed by the 
direction of the AGN jet and its angle with respect to the line of sight. That made it 
appropriate to use the collinear approximation. Here the cosmic rays do not have a fixed 



r & 


8.33 kpc 


Pe 


0.3 GeV cm- 3 


ro 


4 kpc 


a 


0.11 


r s 


35.24 kpc 


Ps 


0.021 GeV cm~ 3 


k e 


10~ 2 GeV" 1 cm" 2 s" 1 sr" 1 


kp 


1 GeV -1 cm~ 2 s _1 sr -1 



TABLE I: The values in this table are from Ref. [45j. 
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FIG. 2: The gamma-ray spectrum from cosmic-ray electrons scattering off of dark matter, for 
various dark matter particle masses M x and masses of the lightest electron-like scalar particle M e . 

incoming direction, thus we want to keep the full kinematics unconstrained, and we integrate 
over the full solid angle. In doing so, we encounter two distinct logarithmic enhancements: 
(i) when the photon is collinear with the final electron, as in |29j , and (ii) when the photon 
is collinear with the initial electron. Notice that this integration is numerically non-trivial, 
as it runs over five independent variables, which can be taken to be the photon angle f2 7 , the 
direction of the final state electron, and the energy of the initial state electron. We adopt 
here a Monte Carlo global adaptive integration strategy, as implemented in the Mathematica 
package. 

We show our numerical results in Fig. [2j where we employ, for the sake of illustration, 
low lightest-neutralino masses (10 and 50 GeV) and small separations between the mass 
of the lightest selectron and the lightest neutralino (0.5 and 1 GeV). The mass difference 
sets the location of the kinematic condition that corresponds to the end of possible resonant 
effects. As a result, we find a sharp cutoff at 0.5 (black line) and at 1 GeV (red and green 
lines). Notice that the drop-off is of about one order of magnitude. Also, notice that the 
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spectrum is significantly harder at energies below the kinematic drop-off, where it is slightly 
harder than E~ 2 . In principle, this might make this particular process distinguishable from 
the diffuse galactic background. Also, from the standpoint of the angular distribution, this 
process is not isotropic, as the corresponding intensity correlates with the product of the 
cosmic-ray density times the dark matter density. Finally, we note that the level of intensity 
we predict for the signal for the models under consideration is only 3-4 orders of magnitude 
smaller than the isotropic diffuse background [IB] . 

While a full-scale simulation of the detectability with Fermi LAT of the signal we cal- 
culated here is beyond the scope of the present analysis, we argue that the level of the 
gamma-ray emission predicted e.g. in the models illustrated in fig. [2] might be potentially 
detectable. At a gamma-ray energy of around 1 GeV, in the same units as in fig. [2j the ex- 
tragalactic diffuse emission is at the level of ~ few x 10~ 4 , and has a very regular power-law 
spectrum [TH]. Considering an angular region of interest of 1 sr, the LAT effective area of 
w 10 4 cm 2 at 1 GeV [16], and an observing time of 3 years, T ~ 10 8 s, this corresponds to a 
total number of isotropic background photons collected at energies around 1 GeV on the or- 
der of 10 5 , and a total number of signal events (for e.g. the red line model of figj2]) of around 
10 2 . This would correspond to a relatively low signal-to-noise (for example, N$ / V Nb ~0.3), 
but in principle an excess corresponding to the sharp spectral feature we predict here could 
be detectable with Fermi LAT. A further handle would be to correlate the detected, albeit 
with low statistical significance, bump with its angular distribution, that should trace the 
product of the dark matter times galactic cosmic-ray density. 



B. Protons 



In the case of protons, the calculation is carried out at the parton level first, with a 
subsequent integral over the parton distribution functions (PDFs) fi(x). We thus consider 
dN - 1 f f f >r-^ d<p d 2 a . N 

dE 7 = r ^ J w x J ^ J dEp J dx ^ x W p d^d% iE -^ (6) 



where the flux of cosmic ray protons 



d<\> , f E p y 



dE p p V GeV 



The calculation is analogous to the one done in the electron case. For numerical purposes, 
it is convenient to make the change of variables E g = xE p , with E g the energy of the 
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FIG. 3: The gamma-ray spectrum from cosmic-ray protons scattering off of dark matter, for various 
dark matter particle masses M x and masses of the lightest quark-like scalar particle M q . 

incoming quark. Doing so, there is no dependence on x in the differential cross section and 
the integrals over x and E q can be done separately. 

As we see in Table [T| the normalization for the cosmic-ay proton flux is two orders 
of magnitude higher than the one for electrons. Naively, one might thus hope to get a 
higher photon signal, but this is not actually the case. The coupling constants enter the 
cross section as a\ R , and because of the fractional electric charge of the quarks, almost two 
orders of magnitude are lost there. This is of course only a simple model-dependent accident. 
The further integration over the PDFs, however, also suppresses the signal, because having 
partons at a given energy requires having much larger proton energies, and the spectrum of 
cosmic-ray protons falls very steeply with energy. As a result, we find that the photon flux 
coming from proton scattering ends up being comparable to the flux from electron scattering 
for the particular case of supersymmetric dark matter. 

Using again for illustrative purposes the case of supersymmetry, we choose three 
sets of lightest neutralino-lightest squark mass pairs in Fig. [3j {M x , M q ) = 
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(10,11), (50,51), (300,310) GeV. Notice that, due to the integration over the parton dis- 
tribution functions, the sharp drop-off we found for the cosmic-ray electron case is now 
smeared out over a broader range of photon energies, but is still noticeable. The overall 
normalization is still in a range which is, for some models, within 3-4 orders of magnitude 
of the isotropic gamma-ray background. 

Notice that we find for the proton case a kinematical condition on the minimum quark 
energy which is completely analogous to the electron case [2S] 

E 

q ~ l-£7 7 /Af x (l-cos0 7 )' 1 ' 

from which we see that the final photon energy cannot exceed M x j (1 —cos # 7 ). That explains, 
for example, the sudden drop to zero at around 5 GeV for the line with M x = 10 GeV in 
Fig.|j 



IV. DISCUSSION AND CONCLUSIONS 



We studied the process of Galactic cosmic-ray electron and proton scattering off of dark 
matter with the emission of a photon in the final state. We argued that the gamma-ray 
flux from this process is, for WIMP models, suppressed compared to gamma rays from dark 
matter pair-annihilation. However, if the latter is in turn suppressed, or nonexistent as 
in asymmetric dark matter models, cosmic-ray-dark matter scattering can offer a unique 
opportunity for indirect dark matter searches with gamma rays. We find that, while for 
dipolar and luminous dark matter the intensity of the signal is well below detectability, for 
WIMP models resonances can enhance the signal to a detectable level in principle. 

We illustrated the expected photon spectrum for the particular case of supersymmetric 
dark matter, and for a few particle dark matter and selectron/squark masses. We found 
that the resonant structure entails a sharp spectral drop-off, at energies corresponding to the 
mass difference between the selectron/squark and the lightest neutralino. The spectrum is 
harder than E~ 2 below the cutoff, and softer, and much more suppressed, at higher energies. 
In conclusion, while this process might be sub-dominant with respect to traditional searches 
for pair-annihilation signatures, it might still offer a unique chance of detecting indirectly 
dark matter with gamma-ray telescopes. 
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